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Abstract The magnetic activity of solar-type stars generally increases with stellar rotation 
rate. The increase, however, saturates for fast rotation. The Babcock-Leighton mechanism 
of stellar dynamos saturates as well when the mean tilt-angle of active regions approaches 
ninety degrees. Saturation of magnetic activity may be a consequence of this property of 
the Babcock-Leighton mechanism. Stellar dynamo models with a tilt-angle proportional 
to the rotation rate are constructed to probe this idea. Two versions of the model - treating 
the tilt-angles globally and using Joy’s law for its latitude dependence - are considered. 
Both models show a saturation of dynamo-generated magnetic flux at high rotation rates. 
The model with latitude-dependent tilt-angles also shows a change in dynamo regime in 
the saturation region. The new regime combines a cyclic dynamo at low latitudes with an 
(almost) steady polar dynamo. 
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1 INTRODUCTION 

Sun-like stars are expected to host hydromagnetic dynamos generating magnetic fields in their external 
convection zones. This expectation is supported by the observed correlation of stellar magnetic activity 
with rotation rate. The physical reason for the correlation is the key role of rotation in dynamos (Parker 
1979). On the Sun, chromospheric Call emission and coronal X-ray emission are correlated with mag¬ 
netic fields. The relative Call and X-ray luminosities are, therefore, accepted as measures of stellar 
magnetic activity. Chromospheric (Noyes et al. 1984; Saar & Brandenburg 1999) and coronal (Vilhu 
1984; Vilhu & Walter 1987; Pizzolato et al. 2003; Wright et al. 2011) activities generally increase with 
stellar rotation rate. 

The increase in X-ray emission, however, saturates for Rossby numbers Ro < 0.1; Ro = P I0 t/T c 
is the ratio of the rotation period to convective turnover time. Coronal activity stops increasing with 
sufficiently rapid rotation (Wright et al. 2011). Saturation in chromospheric activity, though less pro¬ 
nounced, can also be seen in Fig. 6 of Noyes et al. (1984). The observed saturation is an important clue 
for the dynamo theory of stellar activity. It indicates that dynamos do not produce stronger magnetic 
fields when the rotation rate of a star increases beyond a certain limit. A conventional explanation for 
dynamo-saturation is, however, still lacking. 

Literature on stellar dynamos is vast. Brandenburg et al. (1994), Jouve et al. (2010) and Igik et al. 
(2011) constructed dynamo models for stars with different rotational velocities but did not consider the 
saturation problem. Karak et al. (2014) discussed the possibility that a change in the direction of mag¬ 
netic buoyancy from radial to parallel to the rotation axis as the rotation rate increases (Choudhuri & 
Gilman 1987; Schiissler & Solanki 1992; Weber et al. 2011) may cause the dynamo-saturation. They 
concluded that though the change in buoyant rise direction reduces the dynamo’s efficiency, this effect 
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does not suffice to explain the dynamo saturation. More recently, Blackman & Thomas (2015) suggested 
that an increase in differential rotation with angular velocity may cause the saturation. The life-time of 
convective eddies in rapidly rotating stars may be reduced due to distraction of the eddies by rota¬ 
tional shear thus leaving less time for Coriolis force to affect the eddies and reduce dynamo efficiency. 
However, neither observations (Barnes et al. 2005; Collier Cameron 2007) nor modeling (Kitchatinov 
& Olemskoy 2012a) of stellar differential rotation show its considerable increase with rotation rate. 

This paper concerns the possibility that dynamo saturation in rapid rotators is caused by saturation 
in the Babcock-Leighton (BL) mechanism of poloidal field generation. This mechanism was introduced 
by Babcock (1961) and first used in a dynamo model by Leighton (1969). It is related to the observed av¬ 
erage tilt of bipolar groups of sunspots relative to the local parallel of latitude (Hale et al. 1919; Howard 
1996). Due to the tilt, magnetic fields of solar active regions have on average a finite poloidal component, 
which contributes to the global poloidal field upon the sunspot groups decay. There is growing evidence 
for the operation of the BL mechanism on the Sun (Erofeev 2004; Dasi-Espuig et al. 2010; Kitchatinov 
& Olemskoy 2011b; Kitchatinov 2014). BL-type dynamo model of Jiang et al. (2013) reproduces the 
basic correlations observed in solar activity. The BL mechanism can be expected to participate in stellar 
dynamos as well. 

The contribution of the BL mechanism to the global poloidal field formation is proportional to 
sine of the (averaged) tilt-angle. The tilt is believed to result from the Coriolis force and therefore 
can be expected to increase with rotation rate. The BL mechanism obviously saturates when the tilt 
angle approaches 7r/2. The characteristic value of the tilt angle for the Sun is about 5° (Howard 1996). 
Assuming the angle to be proportional to the rotation rate, the tilt angle can reach 90° in a (solar mass) 
star rotating about 20 times faster than the Sun. This is roughly where the saturation in coronal activity 
of solar-type stars is observed (Wright et al. 201 1). 

To assess this idea quantitatively, we use stellar dynamo models, which only differ from a solar 
model by a modification of the a-effect (cf., e.g., Krause & Radler 1980) of poloidal field regeneration. 
Two versions of the a-effect dependence on rotation rate are considered. In one of them, the dependence 
is treated globally by multiplying the a-effect by the factor of ■nu\(ar. J P ^/where a® and Pq are 
the averaged tilt-angle and rotation period of the Sun, respectively, and P rot is the rotation period of a 
star. In the other model, a similar procedure is applied to the Joy’s law of the latitude-dependent tilt- 
angle. The second model shows a wider variety of dynamo behaviours but both show saturation of the 
generated magnetic flux for rapid rotation. The models design is described in Section 2. The results of 
numerical computations for these models are discussed in Section 3. Section 4 summarises our main 
findings and concludes. 


2 DYNAMO MODEL 

The dynamo model of this paper differs little from those used before and, therefore, it is described 
here only briefly. The only difference with Olemskoy & Kitchatinov (2013) is in formulation of the a- 
effect. Other model ingredients like differential rotation, meridional flow and eddy transport coefficients 
remain unchanged. They may also depend on rotation rate but we avoid modifying anything else but 
the a-effect to see the consequence of saturation in the BL mechanism alone which is not disguised by 
other modifications. 


2.1 Dynamo Equations 


We consider the axisymmetric magnetic field B in a spherical shell modeling the external convection 
zone of a star, 


B = e^B + V x 


A 


- sin 9 


( 1 ) 


where standard spherical coordinates (r, 9, (f>) are used, e 0 is the azimuthal unit vector, B is the toroidal 
magnetic field, and A is the poloidal field potential. A similar expression for axisymmetric global flow 



Dynamo Saturation in Rapidly Rotating Stars 


3 


V in a co-rotating frame can be written. 


V = e^r sin9Aflf(r, 9) + p 1 V x 



( 2 ) 


where Ail is the characteristic value of the angular velocity variation within the convection zone, / is 
the normalised differential rotation, p is density, and t p is the stream-function of the meridional flow. 

The dynamo equations are written with normalised variables. Distance is measured in units of stellar 
radius R, time - in the diffusive units of R 2 ^ 1 , where tjq is the characteristic value of the eddy diffu- 
sivity. The toroidal magnetic field is normalised to its ‘equipartition’ strength B 0 for which nonlinear 
feedback on the a-effect becomes essential. The poloidal field potential is normalised to aoB^R 3 ^ 1 , 
where ao is the measure of the alpha-effect. Density is normalised to its value po °n the top boundary. 
The stream-function is normalised to poR 2 V^~ , where Vo is the amplitude of the surface meridional 
flow. The same notations are used for normalised variables as for their dimensional counterparts, except 
for the fractional radius x = r/R. 

The normalised equation for the toroidal field of our dynamo model reads 

8B_ _ V (d£dA_ d£dA\ R m dip 8 f B \ R m dip 8 / B \ 

8t x \9a; 89 89 dx J x 2 p 8x 89 \sin(9 ) x sin 9 89 8x \px J 

V 8 / 1 8{sin9B) \ 1 8 f 8(^rj xB) \ 

x 2 89 \ sin 9 89 J x dx \ ^ dx J ’ 


where 


is the dynamo number and 


a 0 Af IR 3 


Rn 


y 0 R 

90 


(4) 

(5) 


is the magnetic Reynolds number for the meridional flow. The right-hand side of Eq. (3) accounts for the 
toroidal field production by differential rotation, advection by the meridional flow, turbulent diffusion, 
and diamagnetic pumping due to the convective turbulence inhomogeneity (see Kitchatinov & Olemskoy 
2012b, hereafter K012, for more details). The square root of p appears in the last term to allow for 
diamagnetic pumping. 

The poloidal field equation differs by representation of the a-effect but is otherwise identical to the 
solar dynamo model of KOI2, 


8A 

~8t 


xsin0 F(9,P rot 


[ a(x, x')B(x', 9) dx' H- 

J px z sin 9 

Xi 


n_ ■ 0 ^_( 1 dA\ 
x 2 89 \sin0 89 ) 



fd^_8A 

\8x~89 


dip 8A \ 
~89~8x J 


(6) 


The first term on the right-hand side of this equation stands for the non-local a-effect of BL type. 
The kernel function a(x, x r ) in this term differs considerably from zero only if x' is close to the inner 
boundary x-, and simultaneously x is close to the top boundary: 


a(x, x') 
<P b(z') 


4>b (x')(p a (x) 

1 + B‘ 2 {x',9y 

i (1 - erf ((x' - Xb)/h b )), cp a (x) = | (1 + erf ((x - x a )/h a )), 


(7) 


where B 2 in the denominator represents the only nonlinearity of our model, which is the magnetic 
quenching of the a-effect, and erf is the error function. The values x\, = X\ + 2.5/u, and x a = 1 — 2.5 h a 
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ensure smoothness of the d-funotion; h a = 0.02 and h\, = 0.002 in the computations of this paper. 
The a-effect in our model, therefore, describes the poloidal field generation near the surface from the 
toroidal field of the thin near-bottom layer. The factor x sin 9 in the first term on the right-hand side of (6) 
corresponds to the definition of the poloidal field potential A by equation (1). The function F((). P TOt ) 
accounts for the dependence of the o-cffcct on stellar rotation rate. This function will be defined in the 
next subsection. 

The initial value problem for the equations (3) and (6) is solved numerically. The initial field is a 
mixture of dipolar and quadrupolar fields, i.e., the equatorial symmetry is not prescribed. The super¬ 
conductor boundary condition is imposed at the bottom boundary and a pseudo-vacuum condition (zero 
toroidal and radial poloidal fields) at the top. 

2.2 Rotation Rate Dependence of the a-Effect 

The poloidal part of the magnetic flux of a bipolar sunspot group is proportional to the sine of the group 
tilt-angle a. We, therefore, assume the BL type a-effect to be proportional to sin(a) as well. The finite 
tilt in turn is believed to be caused by the Coriolis force action on the rising magnetic loops. The tilt- 
angle is therefore expected to increase with rotation rate of a star. The a-effect changes little when the 
sin(a) passes through a broad maximum near a ~ 90° with an increasing rotation rate. This saturation 
of the a-effect may be the reason for the saturation in stellar magnetic activity. 

The tilt angle in a solar mass star rotating with a period P rot can be estimated as a = a 0 P.j / P rot . 
Two models for the dependence of the a-effect on the stellar rotation rate will be considered. In the first 
model, we use a ‘global’ estimation of the tilt-angle to write 

F{9 , P rot ) = cos 9 sin 2 _ Model I, (8) 

sin Uq 

where a 0 is the mean value of the tilt-angle for the Sun. The sina 0 was included as a denominator in 
this equation in order for the dynamo equations with P-function of Eq. (8) to turn into our former solar 
model (K012) at P rot = P 0 . Then, the value of dynamo number 'D = 3.2 x 10 4 for which reasonable 
results for the Sun were obtained does not have to be changed. In the computations to follow d 0 = 5° 
(Howard 1996). 

In the second model, Joy’s law for the latitude dependence of the tilt-angle (Hale et al. 1919) is 
allowed for. Dasi-Espuig et al. (2010) found a 0 = (0.26 ± 0.05)A and a 0 = (0.28 ± 0.06)A from 
sunspot data of the Mount Wilson and Kodaikanal observatories, respectively, where A is latitude. Ivanov 
(2012) confirmed their results but found a somewhat larger proportionality factor of 0.38 ± 0.03 from 
the Pulkovo observatory database. Stenflo & Kosovichev (2012) infered a still steeper dependence of 
a 0 = (0.56±0.01) sin A from MDI/SOHO magnetograms. They also noted that sin A is a natural choice 
for the fit function in Joy’s law because the Coriolis force, presumably producing the tilts, varies as sin A 
with latitude. The larger tilt-angles from magnetograms are probably explained by the contribution of 
plages (Wang et al. 2015). As the BL type a-effect is more related to sunspots of active regions than to 
smaller magnetic features of magnetograms, we employ Joy’s law in the form a 0 = 0.3 cos 0 to find 

P(0,P rot ) = sin(O.3cos0P 0 /P rot ) - Model II (9) 

for the P-function in Eq. (6). As computations of the solar dynamo with such a profile of the a-effect 
have never been attempted before, we have to define the critical value of T> (4) for the onset of dynamo. 
The computed critical value for dipolar modes in the solar (P ro t = P 0 ) model is T> d = 6.63 x 10 4 (IP = 
7.21 x 10 4 for modes of quadrupolar parity). The solar dynamo is probably only slightly supercritical 
(see discussion in Karak et al. 2015). All computations in Model II were performed with T> = 7 x 10 4 , 
which is about 5% supercritical for the solar case. 

Figure 1 shows latitudinal profiles of the a-effect (9) of Model II for several rotation rates. The pro¬ 
file of cos 9 often assumed in dynamo models is shown for comparison. The a-effect initially increases 
with rotation rate, but then saturates and even changes to a decrease at high latitudes for Prot < P&/ 5. 
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Fig-1 Latitudinal profiles of the a-effect (9) of Model II for different rotation rates. The 
lines are marked by the corresponding values of Pq /P rot . The often assumed cos 0-profile is 
shown for comparison by the dashed line. 


This is the reason for dynamo saturation in subsequent computations. Note the sign reversals in the 
alpha profile for P rot = Pq/ 20. The change of sign is important for interpreting the results of dynamo 
simulations for rapid rotators. 

2.3 Other Model Ingredients 

The differential rotation, meridional flow and diffusivity profile in this paper are the same as in 
Olemskoy & Kitchatinov (2013). 

As mentioned in the Introduction, observations and theoretical modeling both suggest that differen¬ 
tial rotation varies moderately with the rotation rate of a solar-type star of a given mass (Barnes et al. 
2005; Collier Cameron 2007; Kitchatinov & Olemskoy 2012a). The same differential rotation specified 
after the approximation of Belvedere et al. (2000) for the helioseismological rotation law is used for all 
stars. 

The eddy magnetic diffusivity varies little in the bulk of the solar convection zone but drops sharply 
with depth near its base. The diffusivity profile of our model mimics this behaviour: 

r]{x) = r] in + ^(1 - r] in ) ^1 + erf . (10) 

Computations were performed for the following values of the parameters in this equation: r]i n = 10 -4 , 
x v = 0.74 and h v = 0.02. Figure 2 shows the diffusivity profile for these parameter values together 
with the functions (b\, and <fr a of Eq. (7) defining the alpha-effect. The diffusivity drops four orders in 
magnitude near the bottom boundary. A non-uniform numerical grid with very small grid spacing near 
the bottom was used to resolve this sharp variation. 

Single-cell meridional circulation with poleward flow on the top and return equatorward flow near 
the bottom is prescribed (see Fig. 3 and Eq. (13) in Olemskoy & Kitchatinov 2013). The value of R m = 
10 for the Reynolds number (5) is used in the computations, which corresponds to the flow amplitude 
Vq ~ 15 m s _1 and background diffusivity 770 ~ 10 9 m 2 s _1 . With this relatively low Reynolds number. 
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Fig. 2 Profiles of diffusivity and the kernel functions (7) of the non-local a-effect of Eq. (6) 
used in our dynamo models. 


the meridional flow in the bulk of the convection zone is not significant. Only the flow near the bottom 
where diffusion is low (Fig. 2) is important for magnetic field dynamics (Hazra et al. 2014). The same 
value of r/o ~ 10 9 m 2 s -1 was used to convert dimensionless time into physical units. 

As the BL a-effect is related to surface active regions, we illustrate toroidal fields of dynamo models 
by butterfly diagrams for the same near-bottom toroidal flux, 

l 

B(0) = sin0 J (j>b(x)B(x,6) dx, (11) 

Xi 


to which the a-effect of Eqs. (6) and (7) is proportional. Note that the weight function 4 >b of Fig- 2 only 
differs considerably from zero in a thin near-bottom layer. Equation (11) therefore represents the near¬ 
bottom toroidal flux. The factor sin 9 in Eq. (11) accounts for the dependence of the length of toroidal 
flux-tube on latitude. The probability of the production of spots is assumed to be proportional to this 
length. The unsigned total near-bottom flux, 

1 7T 

f m = Jjsin0x(l) h (x)\B( x ,6)\dxd6, (12) 

Xi 0 

is used as a proxy for the overall magnetic activity. 

Dynamo computations started from a weak poloidal field with a mixed equatorial parity. Upon 
several diffusion times, initial field growth saturates. The results of the next section correspond to these 
saturated dynamo-regimes. 
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Fig. 3 Dependence of unsigned flux (12) on rotation period for Model I. The dashed lines 
show maximum and minimum values of the flux in the dynamo-cycles. The full line shows 
the median value (half-sum of maximum and minimum). 


3 RESULTS AND DISCUSSION 

In all runs, the dynamo-fields eventually approached dipolar parity. This seems to be a common feature 
of dynamo models with relatively large diffusivity in the bulk of the convection zone (Chatterjee et al. 
2004; Jiang et al. 2007; Hotta & Yokoyama 2010). The global magnetic field of the Sun is also close to 
dipolar parity (Stenflo 1988). All the simulated dynamos are cyclic. 

3.1 Model I 

Observational statistics of stellar activity usually do not cover epochs comparable to expected activity 
cycles. It is therefore not clear for which phase of the dynamo-cycles the fluxes ( 1 2) should be estimated. 
Figure 3 shows the maximum, minimum and median values of the fluxes in dependence on the rotation 
period. 

The main feature of this plot is that the / m stops increasing with rotation rate for small P rot . 
Observations of stellar activity show its saturation at Rossby numbers 10-20 times smaller than the 
solar value (Fig. 2 of Wright et al. 201 1). Stellar structure in our computations is not varied (because 
the Sun remains the only star for which the tilt-angles are measured). The Rossby number is, there¬ 
fore, proportional to the rotation period. The position of the saturation region in Fig. 3 roughly agrees 
with observations. This figure also shows a decline in / m for still shorter P rot . It is not clear whether 
the decline is real or not because observational statistics have considerable scatter and usually do not 
distinguish spectral types (cf., however, Pizzolato et al. 2003). 

Coronal activity data for moderately rotating stars (Ro > 0.2) are customarily fitted by the power 
law (Pallavicini et al. 1981; Wright et al. 2011), 


Lx oc Ro 7 . 


(13) 
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Fig. 4 Time-latitude diagram of Model I for P rot = 1 day. Top panel: near-bottom toroidal 
flux of Eq. (11). Bottom panel: radial magnetic field on the surface. 


The slope of the plot in Fig. 3 is not constant. However, it varies little in the region of 3 < P rot < 
10 days. The slope <5 defined in the same way as in Karak et al. (2014), 

fl PZt, (14) 

is about S ~ 1.5 in this region. Though our model is designed to qualitatively explain activity saturation 
in rapid rotators and not to achieve quantitative agreement with observations, a preliminary comparison 
is tempting. If 7 is taken to be equal to its canonical value of 2, then our model would mean that Lx 
is proportional to / m to the power of 2.7, which is somewhat larger than the power of 2 expected for 
Lx being proportional to magnetic energy. The power index is still larger (2"f/S ~ 3.6) for the value of 
7 = 2.7 of Wright et al. (2011). 

Measurements of large-scale stellar magnetic fields are most relevant to predictions of dynamo 
models. The measurements are provided by the Zeeman-Doppler imaging (ZDI) technique (Donati & 
Brown 1997). Though statistics of ZDI for main-sequence dwarfs is not vast (about 100 stars), the 
dependence of the large-scale fields on rotation rate has been estimated by Vidotto et al. (2014). Their 
results generally confirm the proportionality of X-ray luminosity to the square of unsigned magnetic flux 
thy of large-scale fields. Vidotto et al. (2014) also found the power law <by oc f?o _119±014 for Ro > 
0.1 and saturation of <1>y at smaller Ro. We have to conclude that our dynamo model underestimates the 
power index 6 of Eq. (14). 

Figure 4 shows the butterfly diagram for the rapid rotator with P rot = 1 day. It is quite similar to 
the solar model (K012). Field patterns computed for different rotation rates are qualitatively similar and 
differ mainly in the field amplitude. This is because only the a-effect in our model depends on rotation 
rate. The cycle period for Fig. 1 is about 10 yrs. The cycle period changes little with rotation rate. This 
is because the meridional flow is not varied. The cycle period is controlled by the flow circulation time 
in advection-dominated dynamos (Charbonneau 2010; Choudhuri 201 1). 
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Fig. 5 Time-latitude diagram of Model II for the solar model. Top and bottom panels have the 
same meaning as in Fig 4 but are computed for the a-effect of Eq. (9) and with P ro t = P©. 


3.2 Modelll 

Model II with the a-effect of Eq. (9) has not been applied to the Sun before. We therefore consider the 
solar case first. Figure 5 shows the time-latitude diagram computed for P rot = P Q . The polar branch of 
the poloidal field is lacking but otherwise the diagram is satisfactory. The cycle period in Fig. 5 is about 
9.6 years. The cycle period varies little with rotation rate (until a new dynamo regime with a not well 
defined period onsets for rapid rotation). The reason is the same as for Model I: the meridional flow is 
kept constant. 

Figure 6 shows the dependence of unsigned flux (12) on rotation rate. The dependence is not as 
smooth as in Model I but generally shows a saturation as well at P mt < 3 days. The power index of 
Eq. (14) S ~ 1.3 (at P rot ~ 7 days) agrees with the computations of Karak et al. (2014). 

The kink in the dependence at rotation periods slightly below 2 days signifies a change in the 
dynamo regime. The new regime is illustrated by the time-latitude diagram of Fig. 7. This figure shows 
field reversals at low latitudes but high-latitude fields do not change sign. Polar fields oscillate about a 
non-zero mean value. 

The following interpretation for this behaviour can be suggested. The positive a-effect (in the north¬ 
ern hemisphere) of solar-type dynamos generates a poloidal field from the toroidal one, and the new 
poloidal field has an opposite sign compared to the old poloidal field from which the toroidal field is 
wound by the (solar-type) differential rotation. Sign reversals of magnetic fields, therefore, take place 
and the dynamo is cyclic. If the a-effect were negative, a new poloidal field of the same sign as the old 
one would be generated and steady dynamos could be expected. Figure I shows sign reversals of the a- 
effect in the northern hemisphere for sufficiently rapid rotation (P rot < 2.2 days). The a-value remains 
positive near the equator but changes to negative at high latitudes. When the polar region with negative 
a is sufficiently broad, a new regime which combines cyclic equatorial with steady polar dynamos sets 
in. 
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Fig. 6 Same as Fig. 3 but for Model II. The two dashed lines and solid line show the maximum 
and minimum values of f m (12) in dynamo-cycles and their median value, respectively. 
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Fig. 7 Same as Fig. 5 but for P rot = 1 day. The fields of high latitudes oscillate about non-zero 
mean values. 
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4 CONCLUSION 

The possibility for the saturation of magnetic activity observed for rapidly rotating solar-type stars 
(Vilhu 1984; Pizzolato et al. 2003; Wright et al. 2011) being related to the properties of the BL mech¬ 
anism for magnetic field generation was discussed. This mechanism naturally saturates when the mean 
tilt-angle approaches 90°. Computations with dynamo models generally support this idea. 

Two dynamo models were considered. The tilt-angle dependence on rotation rate was treated glob¬ 
ally in Model I of Eq. (8). Model II of Eq. (9) employed Joy’s law for latitude dependence of the tilt- 
angles. Both models show saturation of dynamo-generated magnetic flux at high rotation rates. 

Model II is probably more realistic and it demonstrates more variable behaviour. As rotation rate 
increases, the a-effect first saturates, then decreases and can even change sign if the tilt-angle exceeds 
180°. In Model II, this happens earlier at higher latitudes (Fig. 1). This leads to a change in dynamo 
regime with increasing rotation rate. The new regime emerging at high rotation rates combines field 
reversals at low latitudes with oscillations about a non-zero mean value without a change of sign in 
the polar fields (Fig. 7). Magnetic energy is not sensitive to the field sign and, therefore, shows doubly 
periodic variations in this dynamo-regime. The double cycles have indeed been observed in solar-type 
stars with relatively low Rossby numbers (Saar & Brandenburg 1999). 

Two versions of the a-effect are currently discussed for the Sun: the BL mechanism and the a- 
effect of convective turbulence by Parker (1955). Evidence for participation of the BL a-effect in the 
solar dynamo is growing (Erofeev 2004; Dasi-Espuig et al. 2010; Olemskoy et al. 2013). Interpretation 
of dynamo-saturation in rapidly rotating solar-type stars adds one more argument to this line. This is not 
only because this effect saturates for rapid rotation. The a-effect of convective eddies should saturate 
as well when the eddies are twisted by about 90° (the saturation is somehow missed by quasi-linear 
theories of the o-effect). However, the saturation of magnetic activity in the models based on the BL 
mechanism is found at rotation rates about 20 times the solar value where observations also indicate it 
occurs in solar-mass stars. A combined action of both effects is possible (Passos et al. 2014) but the BL 
mechanism seems to be dominant. 

Only the a-effect depends on rotation rate in our model. Other parameters can also depend on P rot 
but were kept constant to see the exclusive consequences of a-effect saturation. Work combining the 
model for (variable) stellar differential rotation and meridional flow (Kitchatinov & Olemskoy 2011a) 
with a dynamo model is currently in progress. 

Acknowledgements This work was supported by the Russian Foundation for Basic Research (project 
no. 13-02-00277). 
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